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Q\ ' We estimate the error in the semiclassical trace formula for the Sinai billiard 

under the assumption that the largest source of error is due to Penumbra diffraction, 
that is diffraction effects for trajectories passing within a distance R ■ 0((kR)~ 2 / 3 ) 
to the disk and trajectories being scattered in very forward directions. Here k is the 
momentum and R the radius of the scatterer. The semiclassical error is estimated 
by perturbing the Berry-Keating formula. The analysis necessitates an asymptotic 
' analysis of very long periodic orbits. This is obtained within an approximation 

originally due to Baladi Eckmann and Ruelle. We find that the average error, for 
sufficiently large value of fci?, will exceed the mean level spacing. 



J??' 1 Introduction 

i 

During the early days of the trace formula nobody really believed that it could be used 
to predict individual eigenvalues, at least not in the strict semiclassical limit. There 
O ■ are mainly two (related) sources of errors. First, the semiclassical energy domain Green 

function is obtained by Laplace transforming the Van-Vleck propagator [jij. However, 
quantum evolution follows classical evolution only for a limited time, a time that seems 
j_j ■ to be longer than first expected but still limited. Computing the Laplace transform 



(with time going from zero to infinity) of such an object is of course adventurous. 

Secondly, the trace formula is obtained by taking the trace of this energy domain 
Green function by stationary phase technique. This is the procedure that selects out the 
periodic orbits. Whether or not this stationary phase approximation is justified for a 
particular cycle depends on h, for sufficiently small fi it is always justified. So, the the 
performance of the trace formula depends on the set of cycles that are required to resolve 
a particular state, and the accuracy of their semiclassical weights. 

Consequently, the semiclassical error depends on the context, the method by which 
the eigenvalues are extracted. In the Berry-Keating method|| it depends only on the 
error of the amplitudes of periodic orbits up to certain length, whereas in complex meth- 
ods, based on cycle expansions the error of the long cycles effects the convergence 
of the cycle expansion and the final position of the eigenvalue. An alternative approach 
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is based on the boundary integral method and the semiclassical limit of the charac- 
teristic determinant ||] . Bogomolny's transfer matrix method is closely related to this 
approach ||]. 

However, it turned out from, numerical computations that the Gutzwiller-Voros zeta 
function does exhibit complex zeroes quite close to the (real) quantum eigenvalues, at 
least in the lower part of the spectrum0, ||, ||]. The Berrry-Keating formula performed 
even better || [l(J. This approach is based on a functional equation for the (exact) spectral 
determinant. By insisting on using the functional equation in the semiclassical limit, one 
actually put in information into the computation: the spectrum is real. The result may 
be a quite respective number of eigenvalues, only a few percent wrong. 

None of these computations indicated what will happen in the strict semiclassical 
limit. Common estimates have suggested that the semiclassical errors, measured relative 
the mean level spacing, should tend to a constant as % —* oo, for system in two degrees 
of freedom. For a nice review, see ref [11]. 



Two common features of chaotic systems may cause problems to the stationary phase 
approximation. The first is intermittency. From a periodic pont of view, it means that 
long orbits does not need to be very unstable, or more precisely, periodic orbit stabilities 
cannot be exponentially bounded with their length. If a long cycle is "relatively stable" it 
means that the corresponding saddle point is not isolated enough from its surroundings. 

Secondly, the complex topology of generic systems make bifurcations abundant (with 
respect to variation of some parameter). So, for a generic systems there are "almost 
existing" and "almost forbidden" orbits all over phase space. For billiards the trace 
integral scans a function sprinkled with discontinuities. Nearly pruned orbits live near 
such discontinuities and the saddle point integrals has to be provided with cutoffs which 
leads to diffraction effects. Even Axiom-A systems such as the Baker map [12| suffers from 
diffraction effects but to a less extent. For smooth potentials one has to face stabilization 
of cycles close to bifurcation. Sufficiently pruned orbits can be included as ghost orbits 
[13 1 but close to bifurcations uniform approximations has to be invoked [14, 15 1. 

There are two pioneering studies indicating that many of the terms included in the 
Berry-Keating sum are way off, the stationary phase approximation behind them is simply 
not justified. Tanner |l(| studied the dynamics close to the bouncing ball orbit in the 
stadium billiard, and Primack et.al. [|l7]] studied the Sinai billiard and orbits scattering in 
very forward direction or sneaking very close to the disk. Both cases deals with systems 
with neutral orbits and as such intermittent. 

The paper [jl7| by Primack et.al. is the main inspiration for the present study. By 
estimating the size of this penumbra and its scaling in energy, they concluded that "the 
semiclassical approximation fails for the majority of the relevant PO 's (periodic orbits) 
in the semiclassical limit". It thus seems likely that the Berry-Keating approach will 
eventually cease to produce individual eigenvalues. Such a conclusion is by no means 
obvious. Simple conclusions can (maybe) be drawn if the system is uniformly hyperbolic. 
But cycles contribute with very different weights in intermittent systems, such as the Sinai 
billiard. If a minority of non diffractive cycles carried a large part of the semiclassical 
weight one could perhaps argue that the trace formula could be saved. But unfortunately 
the situations is the opposite. The cycles being most prone to penumbra diffraction have 
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large semiclassical weights. The goal of this paper is to make this reasoning more precise. 



In a latter study Primack et.al. [11] seems to tone down the importance of Penumbra 
diffraction. They suggest that the semiclassical error in the semiclassical limit is of the 
order of the mean level spacing, or at most diverges logarithmically. It could thus be 
possible to resolve individual states even in the strict semiclassical limit. The study is 
mainly numerical, it is an ambiguous attempt to extrapolate from finite sets of quantum 
states and periodic orbits into the semiclassical limit. Our present study will not support 
their claims. We will indeed find that the error will irrevocably increase beyond the mean 
level spacing. 

It is evident that a numerical study of the semiclassical error with periodic orbits is, 
least to say, difficult. A central tool in our approach will be asymptotic theory for the 



sets of periodic orbits based on an idea of Baladi Eckmann and Ruelle [18|, |19|, |20|, |2l|, |22J. 

What we do is a model study of an intermittent system in two degrees of freedom. 
One should of course be cautious when trying to generalize the result to other system, 
especially to hyperbolic ones. The trace formula is indeed exact for some hyperbolic 
system, such as the Cat map 24] and geodesic flow on surfaces of constant negative 



curvature [25]. The trace formula will probably perform much better in the hyperbolic 
case than in the intermittent. 

The disposition of the article is very much like a cooking recipe. In section |2| we present 
all the ingredients, such as the semiclassical zeta function, the Berry-Keating formula, 
penumbra diffraction, some classical periodic orbit theory and the BER approximation. 
In sec. H we do the actual cooking. We consider the shift of a zero of the Berry-Keating 
formula if a perturbation, due to an error, is added. This simple exercise gives us the 
semiclassical error in terms of finite sums over pseudo orbits. We then relate these 
pseudo orbit sums to various zeta functions. These zeta functions are then calculated in 
the asymptotic theory which we call the BER approximation. In section |I| we present 
the outcome of these exercises. Then follows (section |5|) a round table discussion about 
the validity of the various assumptions and approximations that underlies the results. 

2 Ingredients 

2.1 The semiclassical zeta function 

The starting point will be the Gutzwiller-Voros zeta function ||] whose zeros is to be 
associated with the quantum eigenvalues. It is represented as a product over all primitive 
periodic orbits p of the systems. 

00 / J[Sp/n-n P f]\ 

^)=nn h-W • (1) 

V m=0 \ l iV Pl IV P ) 

where S p is the action along p, [i v the Maslov index and A p is the expanding eigenvalue 
of the Jacobian. To turn this into a Dirichlet series one first expand the inner Euler- 
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product 



A -n(n-l)/2 / A[S p /n-n v l } \ " 

oo 

=nE c p/ InV, "' pf] • ( 3 ) 

p n=0 



where 



A -n(n-l)/2 



C V,n = {'IT - , ' 7- — s (4) 



|A P K 2 n-=i (i-a 3 



If we now expand the product over p, we obtain a cycZe expansion [27] - a sum over all 
pseudo-orbits, that is all distinct combinations of periodic orbits: a = Pi P1 p 2 P2 ■ ■ ■ Vk h ■ ■ ■ 

Z SC (E) = J2C a e^ Sa ~^ , (5) 

a 

where we have defined the quantities 

C« = Y\_ Cp,n p , (6) 
p 

S a = ^ ] n p Sp , (7) 
p 

= n v^v ■ ( 8 ) 

p 

We will restrict ourselves to billiards, the cycle action S p is then given in terms of the 
geometric length S p = L p ■ k where k = y/2E. The units are chosen such that m = h = 1 
and the semiclassical limit ?i — > is replaced by k — > oo. In the following we absorb the 
maslov indices in the amplitudes C a 

Z sc {k) = C a e iLak (9) 

a 

Since we are considering billiards, the redefined amplitudes C a will still be real. 
Note that the size of the amplitudes is 

to leading order, and thus decays fast with n. The zeta function is not seriously affected 
if one restricts the n's to n € {0, 1}. This amounts to retain only the factor m = in 
(||). The resulting type of zeta function is often referred to as a dynamical zeta function. 
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2.2 The Berry-Keating formula 

The spectral determinant for a billiard obeys the functional equation 

D(k) = D(-k) . (11) 
The semiclassical analogue to the spectral determinant is 

D sc (k) = e- inR< - k) Z sc (k) = J2 C a e i{Lak - n * {k)) (12) 

a 

The idea of Berry and Keating j|] was to postulate that this semiclassical determinant 
also satisfies the functional equation (0). This is of course not exactly true Q but by 
insisting on it one can convert eq. (|j~2| ) to a finite sum 



D sc (k) = 2 Yl C " cos(L Q A; - nN(k)) (13) 

a:L a <L BK 

where 

dN ( k ) MAX 

L B k = k—^ • (14) 

N(k) is the mean spectral staircase function. For a billiard it is, to leading order, given 

by 

*r • (15) 

where A is the billiard area. So the cutoff length is given by 

Lbk = -y (16) 

If neutral orbits are present, and they are in the Sinai billiard, their contribution can 

be included in N(k) which is then decorated by oscillation whose amplitude decreases 
with increasing k. 



2.3 The classical zeta function 

Another central object in our investigation will be the (weighted) evolution operator |27||, 
whose action on a phase space distribution function <3?(x) is given by 

Cl$(x) = J w(x,t)6(x - f\y)My)dy . (17) 

The phase space point x is taken by the flow to /*(cc) during time t. w(x,t) is a weight 
associated with a trajectory starting at x and evolved during time t. It is multiplicative 
along the flow, that is w{x,t\ + t^) = w{x, ti)w(f tl (x), t 2 ). If w = 1, the operator just 
describes the classical evolution of the phase space density. 

In e.g. the Sinai billiard the semiclassical zeta function has a branch cut along the negative imaginary 
k-axis and the equation D sc {k) = D so (— k) cannot hold Q 
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In the following we will restrict ourselves to chaotic 2-D billiards, and will use tra- 
versed length L as 'time' variable. 

The trace of this operator can be represented in terms of periodic orbits in two ways. 
First as a sum 

where r is the number of repetitions of primitive orbit p, having period L p , and M p is 
the Jacobian of the Poincare map, its expanding eigenvalue is A p and w p is the weight 
associated with cycle p. 

The trace can also be written in terms of of a zeta function or Fredholm determinant: 

tr C^ = — [ a+l °° e sL ^M- ds . (19) 
The classical zeta function Z w (s) is given by 

p m=0 \ \ P\ p J 

Again we we suffice with the m = factors, and define instead the classical zeta 
function as 

Z w (s) = l[^l- Wp y^j . (21) 

In case the zeta function is entire, the trace can be written as a sum over zeros s 7 of 
the zeta function 

trL L w = Y J ^ L , (22) 

7 

where e SlL can be interpreted as the eigenvalues of the evolution operator. Precisely as 
the semiclassical zeta function, the classical zeta function Z w (s) can be subject to a cycle 
expansion 

4(s)=E a «W e_Si ° • (23) 

a 

We will eventually use the weight w(x, L) to account for diffraction but for the time being 
it is just an arbitrary weight. 

2.4 The BER approximation 



We note that the identity in sec 2.3, holds also after smearing, 



trC L a =YL p f w r J^~ rLp \ = — r i °° e sL ^e^ 2 / 2 ds (24) 

U 4^ V t-i P \det(l-Mp\ 2TTiJ a - ioo Z w {s) 



p r=l 



where 8 a (.) are gaussians of standard deviation a. We will be interested in very gross 
features of periodic sums. Such information is encoded in the behavior of the zeta function 
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Figure 1: Some notations used in the discussion of the one-disk Green function. 



for small s and can be estimated by the Baladi-Eckmann-Ruelle (BER) approximation [ 18 ] . 
The probabilistic approach that underlies this approximation has a long history, see eg 
It was put into the context of Ruelle resonances by Baladi,Eckmann and Ruelle 



in [18 1 . In 22] the formalism was generalized to include zeta functions with 

general thermodynamic weights, and the evaluation of chaotic averages. 

Intermittent systems fluctuate between chaotic and quasi regular behavior. The basic 
idea is to define a surface of section such that all trajectories from the section back to 
itself traverse a chaotic region at least once. The coordinate on this surface of section 
will be denoted x s . 

Let A s (x s ) the length of a trajectory starting at x s and back to the surface of section. 
Further, let w(x s ) be the weight integrated along this segment. We define the (weighted) 
distribution of recurrence time as 

Pw(L) = J w(x s )5(L - A s (x s ))dx s , (25) 

where we have assumed the measure dx s to be normalized: J dx s = 1. The approximate 
zeta function is then given in terms of the Laplace transform of this function 

Z(s) Z{s) = 1 - J Pw {L)e- sL dL . (26) 
2.5 Penumbra diffraction 

A convenient starting point for deriving the semiclassical trace formula for billiards, and 
to study its limitations, is the Boundary Integral Method [29, [L7], H. The eigenvalues of 
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the problem are those for which the following integral equation has a solution 

u(r( S ))=2 ^±(r(s),r(s'))u(r(s'))ds' . (27) 
J s on s 

The integral is performed along the boundary of the billiard. The function u(r(s)) is 
related to the wave function according to the normal derivative 

The boundary r(s) is parameterized by the Birkhoff coordinate s. 
We can write eq (^) symbolically as a matrix equation 

(I-A)U = , (29) 

having a solution when det(I — A) = 0. We can write this determinant as 

det(I - A) = e tr 1 °s( I " A ) = e~ ^ tr ( A ") , (30) 

where 

tr(A") = 

2 n Jd Sl ... ds n ^(r( Sl ), r(s 2 )) . . . ^(rM.'WJ^W^'W) • (31) 

There is considerable liberty of choosing the Green function G(r, r'). In order to 
study the problem of Penumbra diffraction Primack et.al [17] suggested to use the one 



disk Green function (see below) in (|27|). The integral in (|27| ) need then be performed 
only along the square boundary. 

The one disk Green function reads [3C] 

oo 

G( ri ,r 2 ,A0)= l - (Hf\kr l ) + S i {kR)Hf\kr 1 ))Hf\kr 2 )e^ , (32) 

l=-oo 

where fli \z) and H^\z) are Hankel functions and ri, r2 and AO are explained in fig. 
[|. The phase shift function Sf (kR) is defined by 

= -^D^ • (33) 

Using Poisson resummation we get 

oo 

G( ri ,r 2 ,A6)= G^(n,r 2 ,A9) , (34) 

m=— oo 

where 

G( m )(n,r 2 , A(9) = % -J°° (Hf\kn) + S t (kR)Hf ] (k n )) Hf\kr 2 )e il ^ e+2 ™U£ . 

(35) 

The standard semiclassical result is obtained if we 
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1. Retain only the integral approximation G^ m=0 ^ in ( |35l) for the Green function, 
semiclassically this means that classically forbidden orbits such as creeping orbits 
are neglected. 

2. Use the Debye approximations for the Hankel functions. 

3. Compute the integrals by stationary phase. The term tr(A n ) in ( |3l"| ) will now con- 
tain contributions from all periodic orbits with n bounces on the square boundary. 

In fig. H we plot the circle Green function together with its semiclassical limit. Semi- 
classically there is a discontinuity at d = R (where d is the classical impact parameter) 
marking the transition between the lit region and the shadow. In the exact Green func- 
tion there is of course a smooth transition. The interesting thing is that the exact result 
is suppressed as compared to the semiclassical within a distance d cr u, a destructive inter- 
ference between the rays starts already in the lit region and continues into the classical 
shadow. 



This twilight zone was called Penumbra in [17]. In appendix A we show that 

dcrit = (1 + e max (kR))R , (36) 

with 

e max (kR) = c(kR)- 2 / 3 , (37) 

where c depend only weakly (logarithmically) on kR, n, T2 and R. The subscripts max 
was used in the appendix but will be dropped from now on. 

In the penumbra the usually semiclassical look of the Green function is lost. It 
is hard to foresee how what happens when this object is convoluted with itself a la 



eq. (|3l|), in order to get contributions from periodic orbits doing several passages of 



the penumbra. To use the periodic orbit apparatus described in sec we need a 
weight that is multiplicative along the flow. However, we will not use the cycles to 
actually compute the spectrum, we are only interested in estimating the error induced 
by penumbra diffraction, so we can suffice with a rather crude weight. Before we make 
up our mind how this should be achieved, let us discuss the periodic orbits in the Sinai 
billiard. 

The cycles can be coded by associating a (coprime) lattice vector q with each disk- 



to-disk segment, see sec. 3.3. In the limit of small R any such periodic sequence pt = 
qiq2 . . . q n can be realized in the system, except for the rule that two consecutive lattice 
vector may not be identical p0|]. Suppose now that we increase the size R of the disk. 
Some segment of the periodic orbit, say qi, would then eventually need to go through 
the disk which of course is prohibited, the cycle is then said to be pruned. Let's say that 
this happen when R = Ruf- For R < Rbif, there is another cycle p r = qi a qi&q2 • • • qn 
that actually do scatter at the disk when the companion p t just pass by. When R = R^f, 
p t and p r overlap exactly and when R > R^f they are both pruned. This is an analogue 
of the saddle-node bifurcation in smooth potential, and the only source or pruning in the 
Sinai billiard. For a lucid discussion of pruning in billiards, see [|3l]]. 
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Figure 2: The circle Green function versus impact parameter d for fixed values of k = 30, 
R = 1 and r\ = r2 = 3. The full line represents the exact result and the dashed line the 
semiclassical approximation. 



Obviously there is a close connection between pruning and penumbra diffraction. The 
semiclassical weights drops suddenly to zero when the pair is pruned. We have learned 
from our studies of the circle Green function that the quantum pruning is more gradual. 

To estimate the error we will say the pair effectively annihilate each other if they are 
within the transition region discussed above. Not that since |A„ r | 3> |A pt |, we need in 
practice only consider the removal of pt- 

So we define an cycle to be diffractive if it passes the disk within the distance eR (as 
given by (|37|)). The error of the semiclassical weight C p for a diffractive cycle is thus 
defined as 5C P = —C p . A pseudo orbit is diffractive if at least on of the participating 
prime cycles is diffractive. The error of the diffractive pseudo orbit's amplitude is thus 
5C a = —C a 

To exclude diffractive cycles from cycle sums we introduce a multiplicative cycle 
weight w p such that w p = if the orbit is affected by diffraction and w p = 1 if not. 
Exactly, how this is done is discussed in section The associated pseudo cycle weight is 

w a = l[w^ . (38) 
p 

(with the convention that 0° = 1) where n p is restricted to n p < 1. 

One can also consider families of cycles with n passages of the penumbra. For each 
passage the trajectory can either choose to bounce off the disk or not. So such a family 
thus consists of 2™ members and only one of them is unremoved from the calculation 
according to our rule above. This is the one bouncing at every passage of the penumbra, 
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the semiclassical weight of this one is of course very much suppressed and the neglect of 
this orbit is completely neglible. 

The weight now depends on the parameter e. We will skip the index w on the classical 
zeta function (^) and instead denote it Z(s;e). Traces will be denoted tiC^ etc. 

3 Preparation 

3.1 Perturbation of the Berry-Keating zeros 

Let ko be a zero of D{k) as given by ( |13| ) (with superscript omitted): 

D(ko) = 2^C a cos{L a k-iirN{k)) = . (39) 

a 

All pseudo orbit sums are hence forth subject to the cutoff L a < Lbk, which will not be 
explicitly written out in the sums. 

We are interested in how small errors in the amplitudes C a will effect the location of 
this zero. We thus add a small perturbation 

5D(k) = 2j2$C a cos(L a k - mN{k)) , (40) 

a 

where <5C Q is the error of C a , and try to solve 

D(k) + 5D{k) = . (41) 
We then expand k = ko + Sk and consider the solution to 

D'(ko)5k + 5D(k ) +5D'(k )6k = . (42) 
We neglect the last term, it provide higher order corrections and get 

6k ~ D'(k ) ■ (43) 

We will consider the perturbation of a typical zero, sitting at a distance ~ d~ l (where 
d = 4%) from the neighboring zeros. If we assume that the oscillations of D{k) are 
sine- like we can relate the derivative at a typical zero D'{ko) to the mean square of the 
determinant 



D'(k ) = d^2n 2 < D(k ) 2 > . (44) 
We then get for the mean square of the shift 



First we focus on the denominator of eq. (45). We obtain 



<D{k,f>=\2 2 Y J Cl , (46) 
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assuming that cross terms cancel out, cf section |5,3| . 

The average perturbation is, by the same arguments, given by 



<5D 2 (h) >=2Y J 5C 2 a , 



and 



< 



^ 2> i/ 2= ' d-'/V^^F^d- 1 /^ 



Following our reasoning in sec 2.5 we put 

5C 2 a = (l-w a )Cl . 
Recall that w a is either or 1. We can rewrite the function F as 

1/2 

F(k ) = 1 



E a W a C, 



v r 2 ) 

£—ia ^ 'a / 

Recall that we are only considering pseudo orbits such that n p < 1 so 



\a 



a=p 1 p 2 ...p k \ 



IA I" 1 = IA -A A 1-1 



(47) 



(48) 



(49) 



(50) 



(51) 



We can now rewrite the numerator and denominator of ( pOD in terms of the classical 
weights a(e) introduced in section ^ 



a a 

E c 2 = EW°)l 



(52) 
(53) 



The perturbative approach taken in this section is only valid if the predicted value of 
F is small. If F is approaches unity there is no other interpretation than a failure of the 
Berry-Keating formula to resolve individual states. 



3.2 Treating the pseudo orbit sums 

The goal of this section is to relate the pseudo orbit sums (|52| ) and (|53|) to various zeta 
functions. There is an important distinction between these pseudo orbits sums and the 
cycle expansion (f23[) - the occurrence of the absolute values in eqs. (|52|) and (|53|). To 
deal with this we will borrow a trick from ref. p2| . 

Consider the classical zeta function (that is eq. (|2l| ) with higher m factors omitted) 




Consider the cycle expansion of a derived zeta function (note the plus-sign!) 
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All b a are positive, in fact, they are the summands of eqs (|52| ) and ( [52|) 

b a (e) = \a a (e)\ , 

We rewrite this as 



b(L;e)e- sL dL 



where 



b{L-e) = Y^K 5(L-L C 



(56) 
(57) 
(58) 



Conversely, the function b(L; e) is related to the zeta function by means of an inverse 
Laplace transform 

1 



b(L;e) 



Z + (s;e)e sL ds . 



2ni Ja- 

The semiclassical error can now be written in terms of the function b(L; e) 



(59) 



f L ™b(L;e)dL 
V J Lbk b(L;0)dL 



where Lbk is given by eq. (|T|) and e is given eq. (| 
If we now define yet another zeta function 



(60) 



Zo(s: e) 



w 
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we can rewrite Z + (s; e) as 



Z+( 



Z(s-e) 

If the involved zeta functions are entire we can write 

Z 2 (s;e) 

ret' 
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Z(s;e) 



(61) 



(62) 



(63) 



which should be compared with (^). The asymptotic behavior of b(L; e) is (under much 
milder assumptions on the analytic properties) related to the leading zero so(e) 



b(L; e) ~ res 



Z2(s;e) 
Z{s-e) 



oS0 L = Z 2 (so(e);e) L 
Z'(s (e);e) 



(64) 



This leading asymptotic behavior is all we need to evaluate the integrals in eq. (30) since 
we are interested in the asymptotics of the function F. Moreover, as discussed in section 
2.4, the behavior of zeta functions close to the origin is insensitive to fine details in the 
spectrum of periodic orbits, a theory for the large scale structure of periodic orbits exist 
and will be worked out in detail in the next section. 
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3.3 Implementing the BER approximation for the Sinai billiard 

The BER approximation is very well suited for the Sinai billiard, in particular if the 
scatterer is small. The obvious choice of the surface of section is provided by the disk 



itself ||J. 



Consider now the unfolded representation of the Sinai billiard. A trajectory segment 
from the disk to itself can thus be considered as going from one disk, associated with 
lattice vector (0, 0) to some other disk represented by lattice vector q. All segments going 
to q have essentially the same length q = |q| and the following simple expression for the 
distribution of recurrence lengths will suffice for our purposes. 

p(L;e) =Y,a^e)5{L - q) , (65) 
q 

where 

a q (e) = / w(x s )dx s . (66) 
Jn q 

£l q is the set of initial points x s whose target is q. The approximate zeta function is then 

POO 

Z(s;e) = l- p(L;e)e- sL ds = l-J2^(e)e- s ^ . (67) 
Jo „ 



The weight, as defined in sec. |2.5| will be zero, w(x s ) = 0, if the trajectory starting 
at x s and heading for disk q pass some other disk within a distance eR before actually 
hitting q, otherwise it is equal one, w(x s ) = 1. 

If e = 0, then a q is the phase space area corresponding to disk q. 

The computation of p(L; e) will unfortunately be long and boring. It will run in 



parallel with ref [33] in large parts and we will frequently refer to that paper. All error 
estimates in |33| carry over directly, so we will simply omit them below to make things 
a little more transparent. Nevertheless, we will without hesitation display equations as 
strict equalities, but the reader should bear in mind that all expressions are valid in the 
small R limit. 

Consider now a trajectory hitting disk q. The relation between the phase space 
variables and a, see fig ||, and the scattering angle /3 q on disk q is given by 

i? q sin /3 q — R sin a = q sin(<p — q — a) , (68) 

where # q is the polar angle of the lattice vector q, and i? q is the radius of the target disk. 
The argument of sin(0 — # q — a) is small when q is large so we can expand the sine and 
get 

i? q sin/3 q — Rsina = q(4> — #q — ol) . (69) 
It is easy to see that only disks represented by coprime lattice vectors are accessible. 



In ref. [33] we showed that any coprime lattice vector q can uniquely be written on the 
form q = q' + nq c , where n > 2 and q' represents the lattice points closest to the line 
from (0, 0) to q c . We say that disk q lies in the q c corridor, see fig El. 
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Figure 3: Definition of various angles associated with scattering to disk q. 

Actually there are two such neighboring lattice points for each corridor vector q c , one 
with smaller and one with larger polar angle. Below we assume that q' is the one with 
the larger polar angle, the other case is completely analogous, and is accounted for by 
multiplying by a factor of two on some strategic occasions, see below. 

The disk under observation, q, is potentially shadowed only by two disks, namely 
q — q c and q c , see fig. We implement the weight w(x s ) by simply inflating these two 
disks from radius R to (1 + e)R. 

To see how O qc shadows Q q we replace q in (|69|) by q c and put 

/3 qc = tt/2 

i? qc = (1 + e)R . (70) 
O^c = 0q - V(9 9c) 

This gives an equation (in terms of phase space variables a and 4>) for the relevant 
boundary of r2 qc 

(l + e)R-Rsma = q c ({<j)-9 cl + — -a)) . (71) 

9 Qc 

We can now treat ^ q - qc in the same way. We replace q in (|69|) by q — q c and put 

iV qc = (1 + e)R , (72) 

Vqc = #Q + 1/(9 9c) 

and get the equation for the borderline of ^ q - qc 

-(l + e)R-Rsma = {q-q c )({(f>-9 cl - — -a)) . (73) 

9 9c 
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Figure 4: Scattering to disk q = q' + nq c with n = 2 in the corridor q c 



Next we want to change variables from (a, <f>) to (sin a, sin /3 q )), the relation is given 
by eq <M 



By combining equation (|69| ) and (71) we see that the shadowing of f2 qc corresponds 
to the straight line 



q c sin /3 q + (g„ - g c ) sin a = g(l + e) 



1 

R 



(74) 



Similarly, by combining equation (|6^) and (|73|) we see that the shadowing of ^ q - qc is 
given by 

(q ~ 9c)sin/3 q + g c sina = q(l + e) - — . (75) 

If q were not shadowed, f2 q would be given by —1 < sin/3 q < 1 and —1 < sin a < 1. So 
the integral (|66|) is simply the area of the remainder of this square, lying inside the lines 
given by eqs. (|74]) and (ff5|). 

The integration element dx s in (^) over Q q is (to leading order) 



R 

dx s = - — d(sin a)<i(sin/3 q ) 
47rq 



(76) 



It is normalized in such a way that the integral over one octant of the plane is unity. 
So we arrive at the following results 



R_ 

nq 



r i 

(l/2fl-g-eg/2) 2 ■ 
Qc{q-q c ) • 
(l/2R-q c -eq/2) 2 
(q-qc)(q-2q c ) 

I 



l l 

2RT+FJ 
1 



2 It T Q c >l+e/2 

-Am -Qc) 



< q < 
<q< 
<q< 
< q 



l l 

2i? l+e/2 

(2^ + Qc) \+ e /2 

Hm-Qc) 



(77) 



The result can be described in the following way. All disks inside a radius q < 1+772 
are unshadowed and the corresponding trajectory segments not diffractive. Outside this 
horizon the accessible disks are aligned along corridors, each corridor characterized by 
the vector q c , subject to the condition q c < ^ i+e/2 ' Segments longer than | (^ — q c ) 
in a particular corridor q c are always affected by diffraction. Segments longer than |^ 
are always diffractive. 
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Figure 5: The integral of p(L(£)) versus £. The radius is R = 0.1 and e = 0. The 
(full) staircase curve is obtained from eq. ( |65[) w ith amplitudes d q (e) given by flTfl) . The 
(dashed) smooth curve is obtained from eq. (|87,88,39). 



This means that there is a one to one correspondence between corridors (beyond the 
horizon) and accessible disks inside the horizon, so in order to perform the sum (|67|) we 
just need to know the finite number of coprime lattice points inside the horizon, exactly 
how the sum should administrated we be clear below. 

We can now obtain an approximate zeta function by plugging these d g 's into (|67|). 
The integrated recurrence time distribution Jq p(L')dL' = 2q:|q|<L^q i s plotted in figs. 
H and ^. 

However, if R is small, there is a vast number of of coprime lattice points inside the 
horizon, and according to number theory, they tend to be uniformly distributed over the 
plane. As a matter of fact, one can use this uniformity to turn the sum into an integral 
and write down explicit formula for the distribution of recurrence times (|65|). 

Asymptotically, there are 6L 2 /tt coprime lattice points q such that |q| < L in the 
first octant so that the mean density of coprime lattice points is in an asymptotic sense 

1 9 T 

d c {L) = . (78) 



(i) First we consider the case L < ^ T+T72 ■ 

Then d q is a function of q only, cf. eq. (|77|). The distribution function p(L; e) becomes 

R 
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Figure 6: Same as fig. || but with e = 0.2. 



R 12R 
= d c {L)— = —j- . 

We introduce the rescaled length £ = 2RL and rewrite 

, r _, . 12i? . 1 

P (L(0;c) = — 



7!^ 



l + e/2 



(79) 



(80) 



Next we consider the transition region 1/2R < L < 1/R. 
According to eq. ( |77| ) the amplitudes a q depend on the size of corridor q c = |q c | and the 
length of q: a q = a(5 CJ q). The length g is approximately q « q' + nq c = Aq c + nq c , where 
A is a number such that < A < 1, see fig 0. We get 



P(L; e) = 2 ^ ^ a(g c , q)5 a (L - q) 



(81) 



q c n=2 



= 2^ ^ a(q c ,q)S a (L - q c n - Aq c ) , 

q c n=2 

where we inserted the factor 2 to account for both neighbors of q c . 

We will now turn the sums over q c and n into an integral over the density of coprime 
lattice vectors. The parameter A is uniformly distributed in the interval < A < 1 p0| , 
this means that we can just integrate n from n = 2 to oo. 



p(L;e) = 2 [ 
Jo 



l/2R-eL/2 poo 

dq c d c (q c ) a(q c , L) / 5 a (L - q c n) , 

Jn=2 
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Figure 7: The leading zero Zo(e) versus e compared with the asymptotic formula 



rl/2R-eL/2 ]_ 

2 / dq c d c (q c ) a(q c , L)—9 a (L - 2q c ) 

Jo q c 

rmin(l/2R-eL/2,L/2) } 

2 / dq c d c (q c ) d(q c ,L)— , 

Jo q c 



(82) 



where 9(x) is the unit step function. 

Since we are considering the region ^ i+ e /2 < 9 < Ri+e' we nave rnin(l/2i? - 
eL/2, L/2) = L/2. Next we insert the expression for d(q c ,q) from eqs. ( [77|) and <i c (<7c 
from eq. (|78 



(83) 



, T . Z^ 2 , I2q c 1 fl. (1/2R- L-eL/2) 2 . 
p(L;e) = 2 dq c —± -1 - ^ ^-L_ 

Jl+eL/2-l/2R TT q c <?c(L ~ Qc) 



+2 



L-eL/2-i/2R 12qc 1 R (1/2R - q c - eL/2) 5 



lo vr q c irL (L - q c )(l - 2q c ) 

We change integration variable to 77 = 2Rq c and use as before £ = 2i?L. This leaves us 
the following integral to solve 

p(L(0;e) = — f / e/2 ^(1 - (W-^/ 2 ) 2 ) + r^ 2 - 1 d77 (l-77-eC/2) 2 \ 



(84) 



The result of this integral is displayed below (g£ 
Now remains only the case -^j^ < L < ~: 
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Figure 8: The derivative Z'(zq) evaluated at the leading zero versus e. It approached 
unity as e — ► as predicted by eq. (| 



The calculation is completely analogous to the previous case and we get 

rl/2R-eL/2 l 

P{ L {€)'^) = 2 / die d c (q c )— a(q c ,q = L) 

Jo Qc 

_ ^ ,l/2R-*L/2 1 R (l/2R-q c -eL/2) 2 

Jo c ir q c ttL (L - q c )(L - 2q c ) 

_ 24R fi-4/2 drj{l-r]- ei/2) 2 



(85) 



vr 2 Jo 



It is time to summarize the results. It is natural to display the final results in terms 
of the distribution of the rescaled recurrence lengths £ rather than L. We call the distri- 
bution /(£;e) and it is trivially related to p(L; e) according to 

p(L;e)dL = f(£;e)d£ . (86) 

Inside the horizon we already have 

= * e< t^72 (87) 

Beyond the horizon we get, after having performed the integrals (|84|) and (|85| ) 



,,,2 (2e-6e 2 + [4e-3c 2 -26^]ioge+ 

[4£ 2 - 8£ 2 + 4 - 4e£ 2 - e 2 e - 4e£] log(£(l + e/2) - 1) 

+ [4 £ - ^ _ 4 _ 2 ^2 _ £ 2 ? 2 + lQg + e) _ 2 |) 



l+e/2 



< 6 < 2/e 



2 The results of these two integrals can be summarized in one formula, note the absolute value in one 
of the logarithms 



20 



This is valid up to the point £ = 2/e where the function chokes. After that we have 

/(£;e) = 2/e<£ . (89) 

The integral Jq /(£'; e)d£' if these expressions is plotted in figs. [| and ||. We note that 
the statistical treatment of the lattice vectors works surprisingly well, even for such a 
"large" radius as R = 0.1. 

It is also natural to let the zeta functions depend on a rescaled variable z as defined 

by 

s = 2Rz . (90) 

A power series expansion of the the zeta function is related to the moments of the 
distribution /(£; e) 

Z(z; e) = 1 - / e-^f(i- e) = 1 - £ ^f- / £"7(6 ^ , (91) 

m=0 

These moments can be computed from eqs. (|8"7|),(|S8|) and ( p9| ) and are found to be 

1 -e + 0(e 2 loge) m = 

/*£ m f(£-e)#-i ( 1 + °( el °g e )) m = 1 ( 92 ) 
J t. J {!;,€) at < (0(loge)) m = 2 ■ 

^ 0(l/e m - 2 ) m > 3 

The leading zero zq of the zeta function can be computed from eqs. ( pl| ) and (p2|), 
and is found to be 

z = -e + O(e 2 loge) . (93) 
We are also interested in the derivative of the zeta function evaluated at at this zero 

Z , (z ) = l + O(eloge) . (94) 

In figs. and |8| we compare these asymptotic formulas with results from numerical com- 
putation of the BER zeta function as obtained from eqs. (p7|),(p8|) and (|89|). 

Traces will not be a direct concern to us. But, by computing traces in the BER 
framework, we can get an idea of the influence of non leading zeros by computing the 
trace in the BER approximation. That information is relevant for the next section. The 
computation is done numerically by FFT technique. 

The trace formula is in our rescaled units 

1 ra+ioo A 

trC\ = —l e*— log Z w (z)dz . (95) 

2m Ja-ioa dz 

Traces for e = and e = 0.1 are plotted in figs. || and 10. For small £ the trace (in the 
BER approximation) is given by 

trC\ = exp(6£/7r 2 ) - 1 C < • (96) 
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Figure 9: The trace (|95|) for e = (full line). Comparison is made with eq. 
values of £ (dashed line) and eq. (p8|) for large values of £ (dotted line). 
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Figure 10: The trace (p5[) for e = 0.1. Comparison is made with eq. (ptf) for small values 



of £ (dotted line) and eq. (97) for large values of £ (dashed line). 
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If the reader want to verify this, the following hint should be useful: The trace trC\ in 
the range £ < 3 depends only on the behavior of the distribution /(£) in the same range 
(£ < 3) and / is constant in this range. 

In the large £ limit, the asymptotic behavior is given by the leading zero Zq, and thus 

trCf ~ exp(zo£) ~ exp(-e£) , (97) 

see fig |H]. We observe that this asymptotic result settles down very early, that is long 



before the natural scale 2/e, cf. eq. (89) 



For any finite e the zeta function is entire since p(L) has compact support. But when 
e — > zeros will accumulate along the negative real z- axis, building up a branch cut. 
For the limiting case e = this will lead to a power law correction 

trtf~l-^- , (98) 

7T Z £ 

sec fig|. This power law will not be essential in the following. 

We also need to know the value of Z<i{z§\e) as defined in (|6l|). This zeta function 
contain the square of |A p | in the denominator. This means that at z = zq (which is close 
to the origin for small e) the value of the zeta function will be dominated by the shortest 
cycles, which, for small e, will be non diffractive. This implies that ^(zoi 6 ) tends to a 
constant faster than Z'{zq) do, as e — » 0. A simple estimate, based on the methods in 



22] suggests that 

Z 2 (z ,e) = I + 0{R 2 log R) + 0(R 2 )0{e) (99) 



which should be compared with eq. (|9 

4 Results 

We now possess all the tools we need to finally be able to compute the asymptotic limit 



of the error estimate F. To obtain this we fetch from section |3.2 

F 

where 



■ i- f' , (ioo) 



i BK = 2RL BK = AkR , (101) 

cf eq. (|l4|). To begin with we are only interested in the leading asymptotic behavior of 
the function 6(£; e). From section |3.2| we therefore collect 



z [zo, e) 

valid for large values of £. From section [T^ we find 

z e (103) 



23 



Z\zo;e)~l (104) 
Z 2 (z ;e) ~ Z 2 (0,0) (105) 



to leading order in e, for error bounds please go back to section 3.2. Finally, from section 
2.5 we have 

e = c(kR)' 2/3 . (106) 

We are then in the position to compute the large kR limit of the error estimate, which 
we easily evaluate to 

F ~{'-^r) 'I 1 - cmw* ) ■ (107) 

We observe that this function will definitely approach unity which implies that individual 
eigenstates ceases to be resolved. This final collapse will occur on the scale 

kR ~ (cA)- 3 (108) 

Note that if only one symmetry subspace is considered, then A ~ 1/8, so this might 
correspond to a very high energy. 

But how is this asymptotic expression approached. 

Actually we know that F = if kR is less than some critical value which is given by 
£bk = 1/(1 + e/2) (cf. eq. (|96|) ) or written out explicitly 



The solution is given by 



{kR)threshoid = ( l^ 3 - c a- 1 / 3 ) 3 , (110) 



where 

a = ^ + 6y/6(? + 324 /A 2 . (Ill) 

Unstable periodic orbit below this threshold are never diffractive. Of course, neutral 
orbits in this range are subject to diffraction corrections but we do only consider the 
error due to unstable orbits. We are thus unable to estimate the semiclassical error for 
small values of kR. 

To estimate the intermediate behavior we introduce a further approximation of eq. 

i r+ ico z 2 (z\e) , „ r +ico i 



which is highly reasonable, cf. eq. ( |99[) . We replace the classical zeta function Z(z, e) 
by the approximate one Z(z,e), the transforms are again computed by FFT technique 
and the resulting function F(kR) is plotted in fig. 11. In the computation we use 



£ » = U ibU c ' tdz ~ Z2(0 - 0) jL. « ^ (112) 
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Figure 11: The error estimate F(kR) versus kR according to eqs (|112| ) (full line) compared 
with the asymptotic expression (p-00|) (dashed line). A = 1/8 and c is chosen to be c = 1. 



A = 1/8 and an arbitrary value of c = 1. We see is a very steep ascent at the threshold 
discussed above and a fast approach to the asymptotic formula (|107l ). It is likely that 
the individual states ceases to be resolved already here. However, it is hard to make any 
safe prediction regarding the crudeness of our treatment of the penumbra problem, see 
also the discussion of c below. However, the function F does approach unity, sooner or 
later, and it is difficult to find an prevarication of this fact. 



There is an issue which we appear to have forgotten, the constant c in eq. (|127 



is not 
First, 



really a constant, it depends weakly on kR, R, r\ and r2, see fig |l] and appendix] 
this dependence is to weak to be able to alter our general conclusions. Moreover, it is 
not obvious how to implement the dependence on r\ and ri- In the BER application we 
consider disk to disk segment whereas in the study of the circle Green function in sec |2.5| 
we consider (square) boundary to boundary segments. To make the exact connection one 
has to convolute the Circle Green function with itself, the outcome of this operation is not 
obvious. So we used c as if it was a constant when implementing the BER approximation. 
Therefore the "constant" c coming out of the other end of the BER calculation should 
have some residual weak dependence on R and kR. 
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5 Discussion of the validity of the various approximations 
involved 

No chain is stronger that the weakest link. The results presented here is based on a long 
series of approximations and assumptions. Some of them may be readily justified and 
should hardly be controversial but some may seem a lot more crude and one may ask if 
they will allow the chain to break. 



5.1 The correction to the semiclassical weights 

On one hand we argued that penumbra diffraction cannot be accounted for by multiplica- 
tive corrections but on the other hand we needed multiplicative corrections to be able 
to use the machinery of the periodic orbit theory and evolution operators. So we sim- 
ply constructed a multiplicative weight inspired by the penumbra diffraction that should 
be able to provide us with an estimate of the error in the Berry-Keating formula. The 



procedure was discussed at some length in sec. |2T5 . 

One could object that it is too crude to approximate the gradual transition of the 
circle Green function with a step function, and suggest a more smooth weight. This is 
in principle possible, but that would make the calculation in sec. [D^ immensely more 
complicated without changing the result in any significant way. 



5.2 The BER approximation 

It is natural that a theory for the asymptotics of the periodic orbits is hard to check 
numerically. In ref. p(J we compared the exact trace formulas with those of the BER 
approximation for lengths up to roughly the horizon 1/2R for the Sinai billiard. The 
results were encouraging but hardly asymptotic. 

However, a range of asymptotic predictions of the BER approximation appears to be 
correct. It does provide the suggested exact diffusion behavior in the regular Lorentz gas 



(with unbounded horizon) [p4L 22] and the related correlation decays p5[ as well as the 



small radius limit of the Lyapunov exponent [33]. We therefore feel confident that the 
approach works even if more rigorous results are called for. 

We introduced some extra approximations, valid in the small R limit, but our expe- 
rience is that they work excellent even for rather large R. 

One could again raise objection that the diffractive weight is discontinuous and that 



this could cause problems for averages to settle down. But we saw in section 3J3 that 
the effect is just to change the pruning rules slightly, the billiard is as discontinuous 
as before and the fluctuations of the pseudo orbit sums appearing in the numerator 



and denominator (SO) are similar. The integrals in eq. (pCf ) are self averaging and the 



fluctuations of the integrands irrelevant for the estimate of the error. 



5.3 The diagonal approximation 

The diagonal approximation underlying eq. ( |46"1) can be verified assuming that the spec- 
tral statistics is given by Random Matrix Theories fJ6|. The diagonal approximation 
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on the diffractive sum behind ( p7|) is natural, in particular since the majority of pseudo 
orbits are diffractive for small e, but by no means obvious. 

I am grateful to Jon Keating for a very useful discussion. This work was supported 
by the Swedish Natural Science Research Council (NFR) under contract no. F-AA/FU 
06420-313. 



A Stationary phase analysis of the circle Green function 

We will consider scatterings in extreme forward angles so only the second term in the 
integral (^) (m=0) contributes. This means that we only need to evaluate the integral 

G< o >(ri,r2,A0) = % - T S t {kR)Hf\kn)HP {kr 2 )e UAe d£ , (113) 

o J-OO 

where <C A# < ir. Actually, this integral is divergent but as long as we study the 
stationary phase approximation, it serves our purposes. We can still safely use the Debye 
approximation for the Hankel functions Ht{kr\) and Hf(kr2) 



y 7rV z 2 — £ 2 

because kr\ 3> £ and kr 2 3> £■ However, the phase shift function needs a more careful 
analysis. The phase shift function Si{kR) is of unit modulus and we call the phase 

Se{k R) = - H ^ {kR) = e *TW) . (115) 
Hf\kR) 

The Green function now reads 

roc 

G ( -°\r 1 ,r 2 ,A6)= A(r 1 ,r 2 )e l< - iA9 - q '^d£ , (116) 



where the slowly varying amplitude A(ri,r 2 ) is composed of the Hankel functions H^\kr\) 
and H^\kr 2 ). The asymptotics of the integral is determined by the phase function 

<&(£) = _ 7 (A; J R,£)-[ v /(A;ri) 2 - £ 2 -£ &rccos{£ /(kn))]-[^(kr 2 ) 2 - £ 2 -£ arccos(£/{kr 2 ))]+TT/2 

(H7) 

We will now investigate the phase of Si{kR) in detail. To this end we will use the 
uniform approximation for Hankel functions relating the phase of Hankel functions to 
the phase of Airy functions according to |}7]] 

^- Ztm-i ^ ' (118) 
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Figure 12: The function ^i{£) versus £ for fixed values of k = 30, R = 1 and n = r2 = 3 



where 
Defining 

we have for small e 



log( 



t+y/e 2 -(kR) 

kR 



arccos(£/(A;i?)))] 2 / 3 £ < kR 
- ^ 2 - {kRJ 2 } 2 / 3 I > kR 



£ = kR(l + e) , 
-x = (V2kR) 2 / 3 (e + 0{e 2 )) . 



(119) 

(120) 
(121) 



Using this expression for x one gets the so called transition region approximation. 
This approximation fails to yield the Debye approximation as its asymptotic limit but 
the nice thing is that there is a considerable overlap because whenever x <C £ 2 ^ 3 we can 
use the transition region approximation and when i>lwe can use Debye. 

j(x) is a complicated function but one has the following asymptotic expressions [37] 



70*0 



7 (x 



-|(-) 3/2 (l + 0(l/(-x) 3 / 2 )) x 

7T 4 

~ ~2 



,3/2 



-.»"' ~ + 0(x 



-3/2) 



X 



-oo 



+oo 



(122) 
(123) 



so we obtain the Debye approximation when x — * oo as expected. 

The stationary phase condition will simply read ^i{£) = A9 (subscripts denote dif- 



ferentiation). In fig. 12 we plot the function ^e(£) versus £ for some arbitrary chosen 
values of r\/R = r%jR = 3 and kR = 30 together with its Debye approximation. The 
plot gives a picture of how the stationary phase will perform for all AO, just place the 
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ruler horizontally at AO. If it goes below the maximum, it intersects the curve twice. 
The corresponding saddle points corresponds to the direct (larger £) and reflected ray 
(smaller €) respectively. 

We are interested in the location of the maximum of ^t{£). This turns out the lie in 
the region where eqs ( p.2 1| ) and ( |122| ) applies. To leading order in (kR) the location of 
the maximum is the solution to the equation 

xe _| ( _ x)3/2 = 2 - 8 /3 (jU2) -i/3 ( R + R | = ( k Ry l / 3 C( ri ,r 2 ) . (124) 

' - R 2 Jr 2 - R 2 J 




By elementary methods one can show that the solution to this equation lies in the range 



Where in this range the solution lies is completely irrelevant for us. We use this liberty 



and choose e in (|3q ) to be (cf. eq. (|12l[) ) 

-iff 

C(ri,r 2 

where 



e max = 2- 1 /3( fcjR) -2/3 ( log _ c{kR) -Vs _ (126) 



c = c(kR,R,r 1 ,r 2 ) = 2- 1 / 3 I \ log kR + I log 2 - l og (-=JL= 

V 3 3 yvf - i? 2 



+ 



# 2 ,/r 2 - « 2 



2/3 



(127) 
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